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Abstract 

Magnetic monopoles, that are particle-like field configurations with which one can associate a 
topological charge, widely exist in various three dimensional condensate systems. In this paper, 
by making use of Duan 1 s topological current theory, we obtain the charge density of magnetic 
monopoles and their bifurcation theory in a charged two condensate Bose-Einstein system. The 
evolution of magnetic monopoles is studied from the topological properties of a three-dimensional 
vector field. The magnetic monopoles are found generating or annihilating at the limit points and 
encountering, splitting, or merging at the bifurcation points. 
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I. INTRODUCTION 



An elementary particle with a net magnetic charge is an old hypothetical particle called 
magnetic monopole arise in classical electromagnetism and has never been seen in real world. 
Modern interest in magnetic monopole focus on quantum field theory, notably Grand Uni- 
fied Theories and superstring theories, that predict the existence of possibility of magnetic 
monopoles. In 1931, Paul Dirac [l| proposed that the magnetic monopole with an attached 
Dirac string may exist in quantum electrodynamics by their phenomenon of electric charge 
quantization. In 1974, it was shown by 't Hooft and Polyakov [2I that a magnetic 
monopole could be regarded as topological excitations in a quantum field theory due to 
the spontaneous symmetry breaking mechanism. The quantized magnetic charge was inter- 
preted as the topological charge of the magnetic monopole. After 't Hooft and Polyakov's 
works, Duan and Ge |4J studied the rigorous topological expressions of many moving mag- 
netic monopoles, which could not be derived from 't Hooft and Polyakov's theory. It also 
revealed the inner structures of the magnetic charge density current and showed that the 
zero points of Higgs field were point-like source of magnetic monopole. Recently, the theory 
of magnetic monopole has been frequently employed in studying the Grand Unified Theo- 
ries, the phase transitions in the early universe, and the topological excitations in condensed 
matter physics. 

In condensed matter physics, there are also topological objects that imitate magnetic 
monopoles. In chiral superconductor and superfluid, the magnetic monopole excitations 
have been well studied by G. E. Volovik {{J, and such magnetic monopole is the analog of 
Dirac magnetic monopole which combined with two Abrikosov vortices or four half-quantum 
vortices. These vortex lines represent the "conventional" Dirac string. Such Dirac-like 
monopole has been investigated also in ferromagnetic spinor Bose-Einstein condensates jf]]. 
Beside the analog of Dirac magnetic monopole, the " 't Hooft-Polyakov monopole" can also 
be introduced to condensed matter physics js, 3|. In spinor Bose-Einstein antiferromagnets, 
such point-like monopole has recently been worked out by a number of authors 3, |8| . More- 
over, in charged two-condensate Bose-Einstein system, such monopole which has a quantized 
magnetic charge and can be regarded as a real magnetic monopole has been proposed re- 
cently by Jiang [9]. The induced magnetic field of magnetic monopole and their rigorous 
density distribution expression have been deduced by using Duan's topological current the- 
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ory [ljj, [20|. As indicated in above paragraph, magnetic monopole excitations have ahead 
been studied in the context of quantum field theory. Therefore, as pointed out in Ref. 
magnetic monopole excitations in condensed matter offer the exciting opportunity to study 
the properties of magnetic monopoles in detail. Undoubtedly, this will lead to important 
new insights into the general topic of topological excitations in a quantum field theory. 

Furthermore, two-gap superconductivity has drawn great interest recently due to the dis- 
covery of the two-band superconductor with surprisingly high critical temperature MgB2 
[lo| . Two-gap superconductivity is being supported by an increasing number of exper- 
imental reports. Principally, the two gap superconductivity can be investigated in the 
frame of a charged two-condensate Bose system 121 ]. This system is described by 



a Ginzburg-Landau model with two flavors of Cooper pairs. Alternatively, it relates to a 
Gross- P it aevskii functional with two charged condensates of tightly bound fermion pairs, or 
some other charged bosonic fields. Such theoretical models have a wider range of applica- 
tions, including interference between two Bose condensates [3], a multiband superconductor 
14 1, two-component Bose-Einstein condensates [l^], and superconducting gap structure of 
spin-triplet superconductor Sr2Ru04 [16(. Using this theoretical model, two typical topo- 



12 1, and the other 



logical excitations have been presented. One is the knotted vortices [111 , 
is the magnetic monopoles [9|. The main purpose of this paper is to discuss the topological 
properties of the magnetic monopole excitations in charged two-condensate Bose system. 

In Ref. js], Jiang has proposed the magnetic monopole excitation in charged two- 
condensate Bose system, and by using the Duan's topological current theory, the rigorous 
density distribution expression of magnetic monopole has been deduced. The topological 
charges of magnetic monopoles can be expressed in terms of the Hopf indices and Brouwer 
degrees. However, Jiang's conclusions base on a very important condition that the Jacobian 
D(cf)/x) 7^ must be satisfied. When this condition fails, what will happen? In this paper, 
we will investigate the behaviors of the magnetic monopole when this condition fails. 

This paper is arranged as follows. In Sec. II, we give a prime view of the derivation of 
the topological structure of magnetic monopoles. The magnetic monopoles are quantized 
at the topological level and their quantum numbers are determined by the Hopf indices 
and Brouwer degree. In Sec. Ill, we introduce the generation and annihilation of magnetic 
monopoles at the limit point. The bifurcation theory of magnetic monopoles at the first- and 
second- order degenerate points are investigated in Section IV and V, respectively. Section 
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VI is our conclusions. 



II. MAGNETIC MONOPOLE EXCITATIONS IN CHARGED TWO- 
CONDENSATE BOSE-EINSTEIN SYSTEM: A PRIME INTRODUCTION 

In order to make the background of this paper clear, in this section we will give a brief re- 
view of the magnetic monopole excitations in charged two-condensate Bose-Einstein system. 
Firstly, let us consider a Bose-Einstein system with two electromagnetically coupled, oppo- 
sitely charged condensates, which can be described by a two-flavor (denoted by a = 1, 2) 
Ginzburg-Landau or Gross-Pitaevskii (GLGP) functional [12|, whose free energy density is 
given by 

F = + i-A^ x \ 2 + -^-\(hd, - i-A,)m 2 \ 2 + V(* lfl ) + f^, (1) 

2mi 1 c 1 2m 2 c 1 8n 

in which 

v&w) = -b a \^ a \ 2 + °^ a \ A + v m^ 2 + ^t\, (2) 

where r] is a characteristic of interband Josephson coupling strength [17]. The properties of 
the corresponding model with a single charged two-condensate Bose-Einstein system are well 
known. And the relevant field degrees of freedom are the massive coefficient of the single 
complex order parameter and a vector field that gains a mass because of the Meissner-Higgs 
effect. What is very important in the present GLGP model is that the two charged fields are 
not independent but nontrivially coupled through the electromagnetic field, which indicate 
that there should be a nontrivial, hidden topology in this system. However, it cannot be 
recognized obviously in the form of Eq.(TjQ). For working out the topological structure and 



studying it conveniently, we need to reform the GLGP functional. Ref. 
of variables p and xi 2 by 



12 1 introduces a set 



^ a = y/2m a px a , (3) 

where the complex \a = IXaJe* 1 ^ satisfying \xi\ 2 + |X2 1 2 = 1 and p has the following expres- 
sion 

1 l^ll 2 |$9| 2 

p 2 = k— + —)- (4) 

2 mi m 2 



Here p is a massive field which is related to the densities of the Cooper pair. Using the 
variables xi,2 and Pauli matrices a, we define the three dimensional unit vector field n = 
(x,c?x), where (, ) denotes the scalar product and x — (xl Xt)- Then the original GLGP 
free energy density Eq.([T]) can be represented as 

F = ^f(dn) 2 + h 2 (dp) 2 + (-C 2 + V(p,m,n 3 ) 
4 lb 

h 2 c 2 1 

{t[9^C v - d v Cy\ - n ■ d^n x d u n) 2 , (5) 
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where 



8e 

C M = 2ih[xidf,Xi ~ Xldf+Xi ~ X2^%2 + X2 d nX2] ~ —A^. (6) 

Now we find that there exists an exact equivalence between the two-flavor GLGP mole and 
the nonlinear 0(3) a model 18| that is much more important to describe the topological 
structure in high energy physics. In this paper, based on Duan-Ge's decomposable gauge 
potential theory and Duan's topological current theory, we display that there exists another 
kind of topological defect, namely the magnetic monopoles in this system. 

As shown in Eq.([5]), we know that the magnetic field of the system can be divided into 
two parts. One part, the contribution of field C M is introduced by the supercurrent density 
and can only present us with the topological defects named vortices, as what is in the 
single-condensate system. The other part is the contribution n ■ d^n x d v n to the magnetic 
field, which originates from interactions of Cooper pairs of two different flavors and is a 
fundamentally important property of the two-condensate system. 

The induced magnetic field due to n ■ d^n x d u n term is expressed as 

B , = £ ^ £abcn a dvn b dxn c (7) 

87re 

Then, the divergence of the induced magnetic field, namely Q, can be represented in terms 
of the unit vector field n a as 

he 

Q = = e ^ x e abc d^n a d v n b d x n c , (8) 

87re 

this is just the magnetic charge density of the system p m , which is the time component of 
the topological current 

he 

J ™ = 7—^ vXp ta b cd u n a d x n b d p n c } (//, u, A, p = 0, 1, 2, 3). (9) 
one 



It is easy to see that the current is identically conserved, 

3 1^=0 



(10) 



In order to investigate the topological structure of the magnetic charge current, we introduce 
a three-component vector order parameter <fi = <f) 2 , <ft 3 ) formed by the unit vector n which 
satisfies 

n a = ^l U\\ = V^, (a = 1,2,3). (11) 



Obviously, the order parameter can be looked upon as a smooth mapping between the 
three-dimensional space X (with the local coordinate x) and the three-dimensional Euclidean 
space R 3 <fi : x i — > <j){x) G R 3 . n a is a section of sphere bundle S(X). 



Applying Duan's topological current theory [19|, |20|], one can obtain 



Jlk = ^(tiD^), (12) 

and the Jacobian D^{-) is defined as 

e «*wr) = ^ vXp d u <p a d x <p b d p <p c . (13) 
x 

The delta function expression (fT2l of the topological current tells us that it does not 
vanish only at the zero points of <ft, i.e., the sites of the magnetic monopole. The implicit 



function theorem 2l|] shows that under the regular condition 



£°R ^ 0, (14) 

x 

the general solutions of 

<p a (x\x 2 ,x 3 ,t) = 0, (a = 1,2,3). (15) 

The solutions of Eq.( fT5i) can be generally expressed as 

x 1 = x]{t), x 2 = x 2 {t), x 3 = x 3 {t), (i = 1, 2, - • • , K) 

that represent the world lines of K isolated zero points zl(t)(i = 1,2, •• • ,K). These zero 
points are just the magnetic monopole excitations, and the ith world line Zi(t) determines 
the motion of the ith magnetic monopole. 
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The 6— function theory 22j demonstrates the relation 

i=l \ LJ \ x )\zi 

where the positive integer fii is the Hopf index of <fi— mapping, which means that when r 
covers the neighborhood of the zero point ££(t) once, the vector field <fi covers the corre- 
sponding region in space /3j times, which is a topological number of first Chern class and 
relates to the generalized winding number of the <fi— mapping. With the definition of the 
vector Jacobian (j!3p . and using the implicit function theorem, the general velocity of the 
ith magnetic monopole can be introduced 



a z » DM*) 

y/J, _ U4, i _ ^ \ x l 



dt D(Z) 



V° = 1. (16) 



Then, we can get the magnetic charge current in the form of the current and the density 
of a system of K classical point particles in (3 + 1)— dimensional space-time with topological 
charge Wi = (3^ : 

f m = - J^WiV^ir-Ziit)), 

e tr 

Pm = -5\$)D(t) = ^ V W^ir- *(f)), (17) 
e x e 

i=i 

where rji = sgn(D(^)\^ i ) = ±1 is the Brouwer degree, and Wi = (3^ is the winding number 
of 4> at the zero point zl(t). It is clear that Eq. lfTTj) describes the motion of the magnetic 
monopoles in space-time, and the topological quantum numbers are determined by the Hopf 
indices j3i and Brouwer degrees rji of the 0— mapping at its zeros. Here, rji — +1 corresponds 
to a magnetic monopole and rji = — 1 corresponds to an anti-magnetic monopole. 



III. THE GENERATION AND ANNIHILATION OF MAGNETIC MONOPOLES 

As investigated before, the equations of 0's zeros play an important role in describing the 
topological structures of the magnetic monopole in charged two-condensate Bose-Einstein 
system. Now we begin discussing the properties of the zero points, in other words, the 
properties of the solutions of Eq. (Tl5]) . As we knew before, if the Jacobian 

D°A^0, (18) 
7 



we will have the isolated zeros of the vector field <fi. The isolated solutions are called regular 
points. However, when the condition ([18]) fails, the usual implicit function theorem 2lj is of 
no use. The above discussion will change in some way and lead to the branch process. Now, 
we denote one of the zero points as (t*,x*). Let us explore what happen to the magnetic 
monopoles. In Duan^s topological current theory, there are usually two kinds of branch 
points, the limit points and bifurcation points, satisfying 



and 



x 



X 



(t*,x*) 



{t*,x*) 



0. 



1,2,3 



1,2,3, 



(19) 



(20) 



respectively. Here, we consider the case f[T9"j) . The other case (|2T)|) is complicated and will 
be treated in Section 3 and 4. In order to be simple and without lose generality, we choose 
i = 1. 

If the Jacobian 

t^O, (21) 



x 



(t*,x*) 



we can use the Jacobian D l (&) instead of for the purpose of using the implicit function 

theorem. This means we will replace the timelike variable x° = t by x 1 . For seeing this 
point clearly, we rewrite the equations of ( |T5l) as 

<p(x\x 2 ,x 3 ,t) = 0. (22) 

Then we have a unique solution of Eq.( fl5l) in the neighborhood of the limit point (t*,x*) 

t = t(x 1 ), x 2 = x 2 (x l ), x 3 = x 3 (x 1 ), (23) 

with t* = t(x u ). We call the critical points (t*,x*) the limit points. In the present case, we 
know that 

(24) 



dx 1 



(t*,x*) 



i.e., 



dt 
dx 1 



{t*,x*) 



OO, 



(t*,x*) 



0. 



Then the Taylor expansion of t = t(x 1 ) at the limit point (t*, x*) is 



t-t* = - 



2 (dx 1 )' 



[x 



(25) 



(26) 



(t* ,x"* ) 



which is a parabola in x 1 — t plane. From Eq. (l26|) we can obtain two solutions x\(t) and 
x\(t), which give two branch solutions (world lines of magnetic monopoles). If 

dH 



(dx 



1\2 



> 0. 



(27) 



(t*,x*) 



We have the branch solutions for t > t* [see Fig. 1(a)]; otherwise, we have the branch 
solutions for t <t* [see Fig. 1(b)]. These two cases are related to the origin and annihilation 
of magnetic monopoles. 




(a) 




t* * 

(b) 

FIG. 1: Projecting the world lines of magnetic monopoles onto (x 1 — t) plane, (a) The branch 
solutions for Eq. (|26p when d 2 t/(dx 1 ) 2 \( t * t gA > 0, i.e., two magnetic monopoles with opposite charges 
generate at the limit point, i.e., the origin of magnetic monopoles. (b) The branch solutions for 
Eq. (|26p when d 2 t/(dx 1 ) 2 \( t * < 0, i.e., two magnetic monopoles with opposite charges annihilate 
at the limit point. 

One of the results of Eq. (1241) . that the velocity are infinite when they are annihilating, 
agrees with the fact obtained by Bray |23| who has a scaling argument associated with the 
point defects final annihilation which leases to a large velocity tail. From Eq. fT2"4"|) . we also 
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obtain a new result that the velocity field is infinite when they are generating, which is 
gained only from the topology of the vector function <fi. 

Since topological current is identically conserved, the topological charges of these two 
generated or annihilated magnetic monopoles must be opposite at the limit point, i.e., 

Ai% = -PhVia, ( 28 ) 

which shows that (3^ = /3i 2 and 77^ = — r]i 2 , which is important in the charged two-component 
Bose-Einstein system. One can see the fact that the Brouwer degree rj is indefinite at the 
limit points implies and can change discontinuously at limit points along the world lines of 
the magnetic monopoles (from ±1 to =Fl). 



For a limit point it is required that )L„ ^ 7^ 0. As to a bifurcation point 24|, it 



must satisfy a more complex condition. This case will be discussed in the following section. 



IV. BIFURCATION OF MAGNETIC MONOPOLES 



In this section we have the restrictions of Eq.( l20l) at the bifurcation points (t*,x*), 

■0, 



D( 



x 



0, D l 



(t*,x*) 



X 



0. 



1,2,3, 



(29) 



(t*,x*) 



which leads to an important fact that the function relationship between t and x is not unique 
in the neighborhood of the bifurcation point (t*, x*). In our dynamic form of charge current, 
this fact can be seen easily from equation ffT6j) 



dx l 



1,2,3, 



(30) 



(t*,x*) 



which under Eq. f)29p directly shows that the direction of the integral curve of Eq. fl30l) is 
indefinite at (t*, x*), i.e., the velocity field of the magnetic monopoles is indefinite at (if:*, x*). 
That is why the very point (t*, x*) is called a bifurcation point. 

Assume that the bifurcation point (t*,x*) has been found from Eq. (fT5l) and (|29|) . We 
know that, at the bifurcation point, the rank of the Jacobian matrix [|^] is less than 3. We 
suppose 



rank[ 



dx 



1 



(31) 



(t*,af*) 
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and let 



X 



dtp 1 dtf_ 
dx 2 dx 3 

dx 2 dx 3 



(32) 



X* 



which means x* is a first-order degenerate point of the 0— mapping theory. (The case that 
x* is a second-order degenerate point will be detailed in the next section.) From (p 1 = and 
(j) 2 = 0, the implicit function theorem implies that there exists one and only one system of 
function relationships 

x 2 = x 2 (t, x 1 ), x 3 = x 3 (t, x 1 ). (33) 
Substituting fl33l) into (p 1 and <f) 2 , we can obtain 



<p b (t, x\x 2 (t, x l ),x 3 (t, x 1 )) = 0, b = 1, 2 



which give 



5>M = -<^ 



5>M = 
£#4> = - £[20>i + £(#**M 

3=2 3=2 k=2 

3 3 3 

£ $4 = - £[0? O X{ + + £(<4xM] - $0. 



3=2 



3=2 



fc=2 



E*5 

3=2 
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£[2<^Vi + £(#**M 

3=2 



k=2 



where 6=1,2; j, k = 2,3; and 

dx 3 



Xn 



dx 3 
~dt 



5 x 00 



d 2 x 3 
"dt 2 



"> x oi 



d 2 x 3 
dtdx 1 



d 2 x j 



.I'- 



ll 



(<9x 



1\2 ' 



(34) 

(35) 
(36) 
(37) 
(38) 



(39) 



dt 



dx 1 



dx 3 ' 



b _ d 2 <j> b 



9 



2jJ> 



01 dtdx 1 ' 



(40) 



2Ab 







2±b 







2±b 



(41) 



Vu (dx 1 ) 2 ' Vi0 dtdx 3 ' Vjl &rW Vjk dx 3 dx k ' 

From these expressions we can calculate the values of the first and second order partial 
derivatives of (1331) with respect to t and x 1 at the bifurcation point x*. 
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Here we must note that the above discussions do not relate to the last component <p 3 (x, t) 
of the vector order parameter 0. With the aim of finding the different directions of all branch 
curves at the bifurcation point, let us investigate the Taylor expansion of 



F(t, x 1 ) = <p 3 (t, x 1 , x 2 (t, x 1 ), x 3 (t, x 1 )), 



(42) 



in the bifurcation point, which must vanish at the bifurcation point, i.e., 



* 1*^ 



F(t*,x 



0. 



(43) 



From (1421) . the first-order partial derivatives of -F^x 1 ) is 



OF 
~dt ~~dt 



i=2 



<9x J 



3 j dF uif> sr^ 



,J 



dx 1 dx 1 ' ^— ' <9x J 1 
i=2 



(44) 



On the other hand, making use of (132|) . (1351) . (|44|) . and Cramer's rule, it is not difficult to 
prove that the two restrictive conditions in (1291) can be rewritten as 



£>(£) 



(**,«*) 



0. 



(t*,x*) 



(45) 



^( 



(t*,aT*) 



at 



(t*,af*) 



By considering (13"2"j) . the above equations give 



dF 
~dt 



0. 



<9F 



dx 1 



0. 



(t*,ar*) 



(46) 



(47) 



(**,»*) 

The second-order partial derivatives of the function F(t, x 1 ) are easily found to be 

3 3 



d 2 F 



^00 



3=2 



k=2 



(48) 



<9 2 F 
dtdx 1 



i=2 k=2 



(49) 



<9 2 F 



(dx 



fc=2 



(50) 
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which at x* = (£*, x*) are denoted by 



.4 



d 2 F 



dt 2 



B 



d 2 F 



(t*,x*) 



dtdx 1 



C 



d 2 F 



(t*,x*) 



where j, k = 2, 3 and 



3 j3 d 2 (j) 3 
, 0oo " 



d 



a A3 



dt 2 ' vm dtdx 1 ' 



(dx v 



(t*,aT*) 



5 



2^3 



'11 



(dx 1 )' 



(51) 



(52) 



ft 



5 



5 



2j,3 



&rW ' dxWx k ' 



(53) 



According to the Duan's topological current theory, the Taylor expansion of the solution of 
(j) 3 in the neighborhood of the bifurcation point can generally be denoted as 



A(x 1 - z}) 2 + 2B{x 2 - zf)(t -t*) + (t-t 



*\2 



0. 



which is followed by 



and 



dr 1 

A ^ 2 



2B 



dx 1 
~~dT 



c = o, 



C ^ +2B d^ + A = °> 



(54) 



(55) 



(56) 



where A, B, and C are three constants. The solutions of Eq. fl55|) or Eq. (1561) give different 
directions of the branch curves (world lines of the magnetic monopoles) at the bifurcation 
point. There are four kinds of important cases, which will show the physical meanings of 
the bifurcation points. 

Case 1(A ^ 0). For A = 4(B 2 — AC) > 0, we get two different directions of the velocity 
field of magnetic monopoles 



dx 1 



1,2 



-B ± VB 2 - AC 
A ' 



(57) 



which are shown in Fig.2. It is the intersection of two magnetic monopoles, which means 
that two magnetic monopoles meet and then depart from each other at the bifurcation point. 

Case 2(A ^ 0). For A = A(B 2 - AC) = 0, the direction of the velocity field of the 
magnetic monopole is only one 

(58) 



dx 1 



1.2 



A 



13 



t' 



t 



FIG. 2: Projecting the world lines of magnetic monopoles onto (x 1 — t) plane. Two magnetic 
monopoles meet and then depart at the bifurcation point. 

which includes three important situations, (a) One world line resolves into two world lines, 
i.e., one magnetic monopole splits into two magnetic monopoles at the bifurcation point 
[see Fig.3(a)]. (b) Two world lines merge into one magnetic monopole, i.e., two magnetic 
monopoles merge into one magnetic monopole at the bifurcation point [see Fig.3(b)]. (c) 
Two world lines tangentially contact, i.e., two magnetic monopoles tangentially encounter 
at the bifurcation point [see Fig.3(c)]. 

Case 3(A = 0, C ^ 0). For A = A(B 2 - AC) = 0, we have 



dt 
dx 1 



B±^AC _ _2B 



There are two important cases: (a) Three world lines merge into one world line, i.e., three 
magnetic monopoles merge into a magnetic monopole at the bifurcation point [see Fig.4(a)]. 
(b) One world line resolves into three world lines, i.e., a magnetic monopole splits into three 
magnetic monopoles at the bifurcation point [see Fig. 4(b)]. 
Case 4(A = C = 0). Eq. (!55|) and Eq. (l56l) give respectively 

doc n dt _ , . 

— = 0, -— = 0. 60 
dt dx 1 y J 

This case is obvious [see Fig. 5], and similar to Case 3. 

The above solutions reveal the evolution of the magnetic monopoles. Besides the encoun- 
tering of the magnetic monopoles, i.e., two magnetic monopole encounter and then depart 
at the bifurcation point along different branch cures [see Fig. 2 and Fig.3(c)], it also includes 

14 




FIG. 3: (a) One magnetic monopole splits into two magnetic monopoles at the bifurcation point, 
(b) Two magnetic monopoles merge into one magnetic monopole at the bifurcation point, (c) Two 
world line of magnetic monopoles tangentially intersect, i.e., two magnetic monopoles tangentially 
encounter at the bifurcation point. 

splitting and merging of magnetic monopoles. When a multi-charged magnetic monopole 
moves through the bifurcation point, it may split into several magnetic monopoles along dif- 
ferent branch curves [see Fig. 3(a), Fig. 4(b), Fig.5(b)]. On the contrary, magnetic monopoles 
can merge into a magnetic monopole at the bifurcation point [see Fig. 3(b) and Fig. 4(a)]. 
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FIG. 4: Two important cases of Eq. (]59p . (a) Three magnetic monopoles merge into one at 
the bifurcation point, (b) One magnetic monopole splits into three magnetic monopoles at the 
bifurcation point. 

At the same time, the remaining component can be deduced by 

dcc^ a a dec _ j » 

- = 4+xi - J = 2 ,3. (61) 

As previous work, the identical conversation of the topological charge shows the sum of the 
topological charge of these split magnetic monopoles must be equal to that of the original 
magnetic monopoles at the bifurcation point, i.e., 

i f 

for fixed I. Furthermore, from the above studies, we see that the generation, annihilation, 
and bifurcation of magnetic monopoles are not gradually changed, but suddenly changed at 
the critical points. 
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t" 



t 



(a) 



x 




(b) 



FIG. 5: Two world lines intersect normally at the bifurcation point. This case is similar to FIG. 4. 
(a) Three magnetic monopoles merge into one at the bifurcation point, (b) One magnetic monopole 
splits into three magnetic monopoles at the bifurcation point. 

V. THE BIFURCATION OF MAGNETIC MONOPOLE AT A SECOND-ORDER 
DEGENERATE POINT 

In the preceding section we studied the bifurcation of a magnetic monopole at a first- 
order degenerate point. In this section, we investigate the branching process of the magnetic 
charge current at a second-order degenerate point x* = (t*,x*), at which the rank of the 
Jacobian matrix [|^] is 




d<f) 



3-2 



1. 



(63) 



(t*,x*) 



Suppose that 



d<f) 



^0. 



(64) 



dx 3 



(t*,x"*) 
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With the same reasons as in obtaining (133|) . in the neighborhood of x*, from = we 

have the function relationship 



OC OC OC i OC J • 



(65) 



In order to determine the values of the first and second order partial derivatives of x 3 with 
respect to t, x 1 , and x 2 , one can substitute the relationship (1651) into <p 2 (x) = and 3 (x) = 0. 
Then, we get 



Fx(t,x l ,x 2 ) = <p 2 (t,x\x 2 ,x 3 (t,x\x 2 )) = 0, 
F 2 {t,x\x 2 ) = (jf^x 1 ,x 2 ^^(t^ 1 ,x 2 )) = 0. 



(66) 



For calculating the partial derivatives of the function F\ and F 2 with respect to t, x 1 , and 
x 2 , one can take notice of (165]) and use six similar expressions to P7|) . i.e., 



dF c 



dt 



0. 



dF r 



dx 1 



dF r 



(t*,x*) 



dx 2 



c= 1,2. 



(67) 



{t*,x*) 



So the Taylor expansions of Fi(t, x 1 , x 2 ) and F 2 (t, x 1 , x 2 ) can be written in the neighborhood 
of (t*,x*) by 

F c (t,x\x 2 ) « A cl (t-r) 2 + A c2 (t-0(^-^*) 

+ A c3 (t - f )(x 2 - x 2 *) + A cA (x l - x 1 *) 2 + A^x 1 - x u )(x 2 - x 2 *) 



+ A c6 {x 2 -x 2 *) 2 = 



(68) 



where c = 1, 2 and 

1 d 2 F r 



A c i 
A c4 



2 <% 2 
1 d 2 F 



(t*,x*) 



d 2 F 

A c2 = 

dtdx 1 



{t*,x"*) 



A c3 =^ 
dtdx 2 



2 (dx 1 ] 



A 



d 2 F r . 



c5 



(t*,x-*) 



dx 1 dx 2 



A 



(t*,x*) 

1 d 2 F r . 



c6 



(t*,aT*) 



2 (<9x 



2^2 



(69) 



Dividing (I6"5j) by (£— t*) 2 and taking the limit t — > t*, one obtains the two quadratic equations 



of # and dx2 



dt 



dt ' 



~dt 



dx 2 
~dt 



.dx 1 



A cl + A c2 — + A c3 — + A ci (—Y + A c5 



doc doc , , doc > ^ 

~dT~dT + c6[ ~dT } 



(70) 
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and further, eliminating the variable one has the equation of ^ in the form of a 



dx 2 



determinant 



A 14 






A 15 v + A 12 

A u 
A 25 v + A 22 

A u 



dt 



A 16 v 2 + A 13 v + A n 
Ai 5 v + A 12 
A 26 v 2 + A 23 v + A 21 
A 25 v + a 22 



dt 









A 2e v 2 + A 23 v + A 



21 



= o, 



with the variable v 



dxf 
dt 



, which is a four-order equation of 



dx 2 
dt 



, dx \ a , dx \ q / \ o / dx > _ 
ail—) + o 2 (— ) + a 3 (—) 2 + a 4 (— ) + a 5 = 0. 



(71) 



(72) 



rft c?t c?t (it 

Hence, different directions of the branch curves at the second-order degenerate point x* is 
structured. The largest number of different branch curves is four, which means an original 
magnetic monopole with the topological quantum f3rj can split into at most four particles at 
one time with charges Pirji (I = 1,2,3,4) satisfying 



(73) 



VI. CONCLUSIONS 



Our conclusions can be summarized as follows: First, in charged two-component Bose- 
Einstein system, we obtained the dynamic form of magnetic monopole and quantized the 
magnetic charge at the topological level in units of — . The topological quantum numbers 
are determined by the Hopf indices and Brouwer degrees (i.e. the winding numbers), which 
are topological numbers. Second, the evolution of magnetic monopoles is studied from the 

— * 

topological properties of a three-dimensional vector field 0. We find that there exist crucial 
cases of branch processes in the evolution of the magnetic monopoles when D(^) ^ 0, i.e., 
rji is indefinite. This means that the magnetic monopoles generate or annihilate at the limit 
points and encounter, split, or merge at the bifurcation points of the three-dimensional vector 
field 0, which shows that the magnetic monopoles system is unstable at these branch points. 
Third, we show the result that the velocity of magnetic monopole is infinite when they are 
annihilating or generating, which are obtained only from the topological properties of the 

— * 

three-dimensional vector field 0. Forth, we must point out that there exist two restrictions of 
the evolution of magnetic monopoles. One restriction is the conservation of the topological 
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charge of the magnetic monopoles during the branch process [see Eq. ( 128]) and (I62l) ]. the 
other is that the number of different directions of the world lines of magnetic monopoles is 
at most 4 at the bifurcation points [see Eq. (1551) and (I56p ], The first restriction is already 
known, but the second is pointed out here for the first time to our knowledge. We hope that 
it can be verified in the future. Finally, we would like to point out that all the results in this 
paper have been obtained only from the viewpoint of topology without using any particular 
models or hypothesis. 
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